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Abstract 

The important problem of weighted sum rate maximization (WSRM) in a multicellular environment is 



^ \ intrinsically sensitive to channel estimation errors. In this paper, we study ways to maximize the weighted 

sum rate in a linearly precoded multicellular downlink system where the receivers are equipped with a 
, single antenna. With perfect channel information available at the base stations, we first present a novel 

OO ' fast converging algorithm that solves the WSRM problem. Then, the assumption is relaxed to the case 

where the error vectors in the channel estimates are assumed to lie in an uncertainty set formed by 
the intersection of finite ellipsoids. As our main contributions, we present two procedures to solve the 
i intractable nonconvex robust designs based on the worst case principle. The proposed iterative algorithms 

m 

solve the semidefinite programs in each of their steps and provably converge to a locally optimal solution 

of the robust WSRM problem. The proposed approaches are numerically compared against each other 

, to ascertain their robustness towards channel estimation imperfections. The results clearly indicate the 

H ■ 

Cd ■ performance gain compared to the case when channel uncertainties are ignored in the design process. For 

certain scenarios, we also quantify the gap between the proposed approximations and exact solutions. 
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I. Introduction and Motivation 

The capacity limits of even well structured network topologies like, broadcast channels, interfering 
multiple access channels etc. are not yet known [1 1. Some research endeavors have established Shannon 
capacity of few channels where the transceiver nodes are equipped with multiple antennas |2|. Since the 
capacity achieving schemes are mostly not amenable to efficient implementation, several suboptimal alter- 
natives have been proposed. Among them, the foremost example includes linear beamforming techniques 
Q-Q. Nonetheless, such methods have mostly been generally explored under the stringent constraint of 
perfect channel state information (CSI) availability at the nodes of interest. In this study, we devise low 
complexity algorithms for Unearly precoded systems that maximize weighted sum rates in a multicellular 
environment. Such multicellular systems contain both the instances of interference and broadcast channels, 
hence, rendering the problem intractable and NP-hard [6 |. In addition, we relax the constraint of perfect 
CSI availability at the base stations (BSs) and propose computationally efficient algorithms that take care 
of the unwanted and unavoidable channel uncertainties. We focus on the case where channel errors are 
contained in a set formed by an intersection of ellipsoids. The channel uncertainty model considered in 
this paper is general enough to cover all models known in the literature. 

Consequences of Ignoring Channel Uncertainty: The deleterious effects of channel errors have been 
noted in earlier studies on classic CDMA systems Q, lUl, where, for instance, using a so-called system 
sensitivity function a capacity loss of as much as 97% has been reported owing to system imperfections. 
For a more contemporary outlook, we focus on the much explored codebook based limited feedback 
systems ll9l- |[TTl . We will base our discussion mostly on the type of limited feedback schemes considered 
in 191, lITOl . and in our case we focus on a linearly precoded point-to-point multiple-input single-output 
(MISO) transceiver. For a given estimate of the channel, h, at the receiver, the index of the codebook 
vector, w*, that maximizes the received signal-to-noise ratio (SNR) is determined. It, along with the 
request of a particular modulation and coding, is then fed back error-free to the transmitter where the 
beamformer corresponding to that index is employed to transmit the data modulated as per the request. 
We now assume that the channel is not perfectly estimated at the receiver and it is corrupted with some 
error, i.e., h = h + 5 where h is the true channel at the receiver and 5 ~ CA/'(0,I) represents the 
errorsj^ With this assumption, we note that the probability of not exceeding calculated SNRopt, where 
SNRopt represents |hw*p, is straightforwardly approximated by l — exi){—^e{X^p^SNRopt — c)), where 

'We associate the Gaussian probability law for the error term for the sake of demonstration and further, assume zero cross- 
interference and unit variance noise at the receiver. 



c is a constant dependent on the known channel estimate, /i^ is the mean of the exponential variable 
corresponding to the first component of the error vector and Xopt represents the only eigenvalue of 
W* = w*w*^. It is interesting to observe that this approximation is close to or exactly (the desired 
event) only when SNRopt ~ Aoptc or c > X^p^SNRopt- For the other scenario c < X^p^SNRopt, the 
probability of exceeding the threshold SNRopt is significantly smallj^ Hence, in the presence of channel 
uncertainties, even with optimally designed codebooks, the above crude calculations show that for certain 
events of practical significance we have a relatively small probability of exceeding SNRopt- This in turn 
implies that the initial request of the receiver for a particular modulation format will not match with the 
actual requirement, and, hence, owing to this mismatch there could potentially be a drastic increase in 
the probability of making decoding errors. 

As a possible remedy for the curse of channel estimation errors, the receiver should take into account 
the true channel. There are two possibilities to model this problem. One includes associating a probability 
distribution (or more generally a family of probability distributions) with the error terms and then 
translating the problem into having stochastic constraints. The stochastic version of the problem is often 
difficult, if not impossible, to solve exactly |12|. Moreover, it requires the knowledge of distribution of 
error terms which, in most practical cases, is often unknown. As an alternative, the second possibility is 
to assume that irrespective of the probability law the error terms follow, they lie in a certain bounded 
region defined by an uncertainty set^ The problem can then be modeled to satisfy the constraints for all 
error realizations. This gives rise to the philosophy of the so-called robust optimization^ Contrary to the 
stochastic version, the worst case robust version of a problem is often tractable or can be approximated 
by tractable set(s) of constraints llT2l . Coming back to our limited feedback example, by adopting the 
robust optimization principle, we may maximize the performance metric (SNR in the present case) over 
all true channel realizations to determine the optimal index at the receiver. Once such a problem is solved, 
it is easy to see that SNRopt can be exceeded in practice for a large fraction of errors. 

Despite the huge interest in the WSRM problem, it mostly remains unsolved in typical scenarios of 
interest. For instance, a linearly precoded system in a multicellular environment that achieves capacity in 
the downlink is yet to be characterized. The main 'culprits' in accomplishing this goal are the broadcast 
and interference channel components that constitute the whole system. For both of these channels, the 

^We numerically quantify a similar negative impact of channel errors via distribution functions in Sec. |V] Fig. |6] 
^ For tractability reasons the set is often assumed to be convex compact. 
''We also follow this strategy for the problem considered in this paper. 



capacity limits are not known although some progress has been made in recent research endeavors | ]T3l . 
|[T4l . In these and related references, the problem setup generally caters for very specific cases either in 
terms of the network topology or in terms of the assumptions made on the signals and systems involved. 
In a similar way, suboptimal solutions for WSRM problem with linearly precoded transmitters have been 
presented for perfect CSI availability in |l3], lITS I and the references quoted therein. A successful recent 
attempt towards characterizing the capacity region in a multicellular environment includes replacing the 
sum rate functions with the surrogate of degrees-of-freedom (d.o.f.) in the so-called interference alignment 
strategies |[T6l . fTJ] . Nonetheless, it is well known that there could potentially be a substantial gap between 
the exact capacity at finite SNR and the d.o.f. achieved for a certain network setup. On top of aU this, 
most of studies conducted to explore the capacity and/or achievable rate region for linearly precoded 
systems, for instance, have nearly always assumed perfect channel estimation. 

In this work, we consider the WSRM problem in the downlink of a multicellular setting with linear 
precoding by relaxing the stringent assumption of perfect CSI. In particular, our contributions include: 
a) a novel polynomial time algorithm that approximately solves the NP-hard problem of WSRM under 
the scenario of perfect CSI; b) the study of the case of CSI corrupted with errors in an affine manner, 
when the errors are present in an uncertainty region formed by an intersection of a finite number of 
ellipsoids; c) a first robust solution to the WSRM problem based on the worst case principle of design, 
where in we approximate the exact robust counterpart following a type of Lagrange relaxation so that 
any solution feasible for it is also feasible for the exact robust counterpart; d) a second solution, where 
we first approximate the uncertainty set of intersecting ellipsoids with an inner ellipsoid of maximal 
volume, and are able to transform the robust counterpart into a tractable equivalent formulation by using 
the well-known 5-lemma |[T2l . ifTSi ; and e) numerical results where we compare the approximations with 
exact solutions in specific scenarios and demonstrate gains associated with both robust solutions to the 
WSRM problem compared with the nonrobust version. 

Organization: Section JI] formulates the problem and details the assumptions made. The solution for 
the case of perfect CSI is presented in Section Hill Sections |IV] delineates details of two approaches to 
solve the WSRM problem with channel uncertainties. Numerical results and conclusions are presented 
in Sec. |V] and Sec. |Vll respectively. 

Notation: Boldface uppercase (lowercase) letters are used for matrices (vectors). The notations for real 
and complex matrix (vector) spaces are conventional. The size (dimension) of vectors and matrices are 
mostly inferable from context or is explicitly mentioned. All unconventional notations are defined in the 
text at their point of appearance. 



II. Problem formulation 

Consider a system of B coordinated BSs and K users. Each BS is equipped with T transmit antennas 
and each user with a single receive antenna. We assume that the data for the kth user is only transmitted 
by one BS. To keep the representation general enough we use set notation to represent the users that are 
served by a BS. Such a notation can, for instance, cover the scenario when a BS is to schedule users 
based on some priority. The set of all users served by BS b is denoted by Ui,. We further assume that 
the cai^dinality of the set Ui, is Kh i.e., Ki, = \L(b\ for all b so that J2b-^b ~ ^- "^^^ tuple (6, k) provides 
the index of the A;th user being served by the 6th BS, and we use B = {1, . . . , B}. Under frequency flat 
fading channel conditions, the signal received by the kth user served by BS b is 



yb,k 



^bi,k^b,kdb,k + X ^bi,k^i,pdi,p + nb,k (1) 



{i,pmb,k) 



where h^. ^ e C^^"^ is the channel (row) vector from BS i to user k being served by BS b, 'Wb,k £ C^^^ is 
the beamforming vector (beamformer) from BS b to user k, db± is the normalized complex data symbol, 
and rifo fe ~ CA/'(0, cj^) is zero mean circularly symmetric complex Gaussian noise with variance u^. The 
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total power transmitted by BS b is ZlfcewJI^^.'^L- "^^^ SINR ^b,k of user k is 

_ \"-bb,k^b,k\ _ \"-bb,k^b,k\ 

^b,k - — - 72 - B '^^^ 

' \"-bi,k^i,p\ 2 v-^ I |2 ll |2 

ig|;^W. ^ + + 2^ }_^\hb^,k^n,l\ 

{i,p)^{b,k) jeUb\k n=l,n=ibleU„ 

where the interference in the denominator is divided into intra- and inter-cell interference power terms. 
We are interested in the problem of weighted sum rate maximization (WSRM) under a per-BS power 
constraint, which for the case of perfect CSI is formulated as 

maximize J2beB J2keU, <^b,k log^il + 7b,k) (3) 

where Ob^k S ^++- The WSRM problem is a challenging nonconvex problem. In fact, recently it has 
been shown to be strongly NP-hard even in the case of perfect CSI availability 0. In order to arrive at 
a tractable approximation of the above problem when the channel information has been corrupted with 
errors, we will first develop a novel approximating algorithm of this program in the case of perfect CSI. 
In the next stage the proposed approach will be further leveraged to the case of imperfect CSI. 



A. Uncertainty Modeling 

The errors in the traditional channel estimation processes are known to follow Gaussian distribution 
|[T9]| . This, of course, is a consequence of the assumption of ignoring other impairments incurred in the 



process. Further, we recall that the most significant probability content is concentrated around the mean 
of a standard Gaussian model and '3-0"' rule is a well accepted manifestation of this factl^ Leveraging 
the same theme to higher dimensional representation of Gaussian distribution, we observe, a similar 
argument reveals that focusing on Kexp(— xRx^) > r should suffice for all practical purposes for a 
properly chosen rl^ This clearly motivates for an ellipsoidal uncertainty set, which in addition to some 
theoretical justification also offers various computational benefits on account of its convexity. 

Traditionally, the uncertainty has been assumed to lie in a given ellipsoid lHOl, IHbQ Nonetheless, such 
an assumption limits the possibilities of modeUng uncertainty in several other cases of practical interest. 
For example, modeling the uncertainty set as a polyhedron is better (compared to a single ellipsoid) 
when the errors are predominantly due to quantization effects lIlTl . Herein, we consider a more general 
case when the uncertainty set is the intersection of elhpsoids. Specifically, the uncertainty in the channel 
vector hfo^ is modeled as 

h6„,fc = hb,^,A; + 6fe„,A,., 5f,„,fe e = {5f,„,fe : 6fe„,fcPfe„,fc5?„,A; < Pb^,k, 0' = 1, 2, . . . , Q} (4) 

where h^^ ^ is the estimated (known) channel and the uncertainty set Sb„,k is composed of an intersection 
of Q ellipsoids^ The above model can be used to mathematically represent different types of uncertainty 
sets. Some examples may include: i) when 65^^^?^ ^^.6^ ^ = [^b„,k]q^b k ^^'^ dimension of the 
vector Sb^^k is Q, we have an uncertainty box {| [5fe„,A,]g| < ^J{(^b k^'^Pbrt^k^Q} representing the error 
region, where we recall that the notation [6b„,fc]g gives the qth component of the vector 6;,^ j^.; ii) when 
the matrix k ~ ^1 k^l J^^*- '•^^ outer product of the column vector fc> we have the polyhedral 
uncertainty set {|6fe^ < y^p^T^, Vg}; and iii) when g = 1 we have the conventional single ellipsoid 

error model. We note that the uncertainty in the channel h^^ ^ can be modeled on similar lines. 

B. Optimization Problem Modeling 

To mathematically formalize the robust principle of Sec. U we consider a function /e(x, z), where 
X G A" C C" is the decision variable and z is the data parameter. For the sake of argument we assume 

^Here denotes the variance of a standard normal distribution. 

*K denotes an appropriate constant tlrat ensures unit area under the multidimensional Gaussian probability density function. 
X is assumed to be a row vector here. 

^We recall that an ellipsoid is a convex set f (Q,c) = {6 : (6 — c)Q(6 — c)^ < p}, centered at c and parameterized by its 
radius p and a positive definite orientation matrix Q. 

^We make some additional assumptions on the uncertainty set as outlined in the derivation provided in the Appendix. 



that the data parameter is perturbed and z G ^, where Z is some tractable uncertainty set. We are deaUng 
with the problem of maxxeA" /e(x, z). As outlined in the introduction, we would like to ensure that the 
function is maximized over all instances of z G Z, or maXxgA' minzg^ /e(x, z). This model of the 
robust, uncertainty immune, optimization problem is dubbed worst case robust counterpart of the original 
problem and this strategy will be adopted in the discussion to follow when we deal with an optimization 
problem affected by uncertainty. We note that this policy has been introduced and popularized recently, 
see |[T2]| and the references therein. Indeed, the proposed approach can be quite conservative, thereby 
leading to pretty diminished objective value. Nonetheless, the philosophy has the additional advantage 
of being unaware of the statistics of the uncertainty vector. 

III. Solution for Perfect CSI 

As mentioned above the optimization program (|3} in its original form is nonconvex and NP-hard. 
Further, it does not appear possible to find an equivalent convex formulation of the problem by, say, 
some substitutions etc. Hence, we need to find the approximate solution of the problem. For this purpose 
we first note that owing to the monotonicity of the logarithmic function, (O can be equivalently cast as 



maximize 



n n ib,k (5a) 

beB k&Ab 

subject to -—- ^ ^ > {t,[^ ■ - 1), 

yb£B,kel(b, (5b) 

\\^b,k\\l < Pb, yb, (5c) 

keUb 

The above formulation is still not amenable to providing us with a solution to the original problem. 
Therefore, we proceed further and again obtain the following equivalent formulation of the above problem 

maximize 1 | 1 f tft a; (6a) 

vrb,k,tb.k,fJ.b.k ' 
b€B keUb 



subject to hb,,k^b,k > y {t^j^k"'" - 'i^)fJ-b,k, lm.{hb,,kWb,k) =0, Vbe B,k e Ub, (6b) 

B 

cr^ + ^ \\lb,,k'^b,j\ + ^ ^ |hb„,A:Wrt,«| < /Ufo^fc, ^b e B,k eUb, (6c) 

j&At\k n=l,n^bl&A^ 

\\^b,k\\l < Pb, V6, (6d) 

In the above formulation we note that the constraint, lm(h.b^,k'^b.k) = 0, is without loss of generality. 
It is due to the fact that a phase rotation of the beamformers does not effect the objective of the 



problem. Similar arguments have also been used in ||221 . Next we note that the constraint in (l6c]) is 
SOC representable. Indeed, = {ij,i,,k — 1)^ — (/^fe.fc + 1)^ implies 

2 ^ 2 1 1 

which is an SOC constraint. Now, we deal with the only nonconvex constraint in the inequality of (|6b] l. 
The troublesome bit in this constraint is the nonconvexity of the function on the right side in the variables 
involved. It is seen that the solution of the optimization problem in Q is invariant to any scaling in fc. 
Thus we can consider the case when ab.k is greater than 1 for all b, k. With this assumption, the function 
^b^k'" '' becomes concave. It is a well known result that the geometric mean of nonnegative concave 
functions is also concave |23|. Therefore, being a geometric mean of ~ 1) ^^'^ f^b,k^ the right 



side of (I6bl ) is a concave function of the two variables. 

Before proceeding forward with the solution, we note that the iterative nature of the proposed approach 
is similar in spirit to the recent work in [15 1. However, the way auxiliary variables have been introduced 
renders it novelty with respect to the earlier work in [15 |. To deal with the nonconvexity of (l6bl ) we will 
resort to a recently introduced sequential approximation strategy in li24l . Summarizing briefly, for each 
iteration the philosophy involves approximating the nonconvex function with a convex upper bound of an 
auxiliary variable such that the gradients of the original function and the approximation are equal for a 
properly chosen additional variable. Mathematically, let F(x) be the function that induces nonconvexity. 
For the kth step, the technique of 12411 involves determining a convex upper bound -Fc(x, y) of the function 
F(x) such that for an appropriate y = /(x), the following relations hold 

F(x) = F,(x,y), VF(x) = VF,(x,y). (8) 

Under the conditions mentioned above, a natural choice for the value of y in the k + 1st iteration is 
Yk+i = /(xfc). Fortunately, being a concave function, the appropriate upper bound of the geometric mean 
on the right side of (|6bl ) is just a first order Taylor expansion i.e.. 



b,k ^ ( {n)s . 

— {^^b,k - 1^1,1)+ 



\ t^b,k 



'^Oib^k 



'b,k 



(n) 

m.k-t'-:i)^f(-\i,f,) (9) 



%,k 

where the superscript n on the right side is used to indicate the value of the approximation in the 
nth iteration of the algorithm to be outlined later. In addition, it is easy to see that the update in the 



n + 1st iteration follows the straightforward rule (^IT^^ '° ' /^il^^^) ~ (^i'fc ' /^i'fc )■ Clearly, the 
conditions mentioned in ([8]l are satisfied for this update function. With this, we have seen the way the 
problem in © can be transformed into a convex optimization framework. We only need to deal with 
the objective in (l6ab . Although not immediately obvious, it is also expressible as a system of SOC 
constraints. For this purpose we recall the result that a hyperbolic constraint of the form < xy is 
expressible as \\{2z, {x — y))^||2 < {x + y), where x, y G M^. Now, by collecting a couple of variables a 
time and introducing an additional squared variable, we can use the SOCP representation of the hyperbolic 
constraint and end up having J2b^b three dimensional SOCPs |[T5l . |[25l . 

Hence, now we can present a convex formulation of the WSRM problem when perfect CSI is available. 
The problem in (|6]l can be approximated in the nth iteration as 



maximize 




(10a) 



subject to hb^^ky^b,k > f^'^\i,il), lm{h.b^,k^b,k) =0, \/b e B,k £ Kb, (10b) 
T^ + Y^ \hb,,ky^b,jf + J]]|hfe„,fcW„,z|^+-(;Ub,fc + 1)2 < -{fj,b,k - lf,yb£B,k£Ub, (10c) 

\WbAl < Pb, V6 (lod) 

keUb 

where the notation (•)soc indicates that the objective admits SOC representation. The algorithm outlining 
the evaluation of the above problem is sketched below 

Initialization: set n := and randomly generate (tlk^rbD- 
repeat 

• Solve the optimization problem in ([TOl l and denote the optimal values of 

{tb,k,f^b,k) by (tb,fc,Mfe,fc)- 

• Set (^iJ'^^/^iT^^) = (.tb,k^ ^^b,k) and update n := n + 1. 
until convergence 

(n) 

It is significant to note that t j^ ^ may tend to 1 for some n, and, thus, may induce a singularity in the 
right side of the inequality constraints given in (llObl i. Hence, it is practical to implicity assume that 
tb > 1 + £, where e > or to regenerate an independent sequence of t j^ ^ on the emergence of such an 
event till an appropriate sequence is found. The convergence proof of the algorithm is given in Sec. IIV-CI 



IV. Robust Designs for WSRM Problem 

Based on the worst case strategy for robust optimization outlined above in Sec. III-B[ determining the 
exact robust counterpart of the WSRM problem given in <^ appears an intractable problem. Hence, we 
will resort to determining the robust version of the approximate solution for WSRM problem developed 
in Sec. |lIIJ To be specific, the exact robust counterpart of the formulation derived previously is given as 

maximize 1 I 1 I (Ha) 

beB kGUb 

subject to |(hb,,fc + 5;„,A:)wf,,fc| > /(i, A), \/b £ B,k £Ub,y5b,,k & <Sb,,k (lib) 

B 

\ihb^± + h„k)wb,jf + ^ '^\0^b„,k + ^b„,k)^n,i\'^< fJ'b,k, yb £ B,k £Ub, 

jeUb\k n=l,n^bl&U„ 

VSfo^^fc G 'Sb,,k,^b„,k G Sb^^k (lie) 

5] l|wb,fc||i < Pfc, V6 (lid) 

keUb 

where the uncertainty sets Sb^,k represent an intersection of ellipsoids and h;,^^ ^ denote the known values 
of channels. The above formulation, although the robust counterpart of the approximate solution in (ITOl l is, 
indeed, still challenging. We observe that in addition to being nonconvex, it also suffers from tractability 
issues. In particular, a part from the power constraints, all remaining constraints are semi-infinite in nature. 
It does not appear possible to arrive at an equivalent tractable version of the above optimization problem. 
Therefore, we need to derive approximate solutions. In what follows, we will present two approximation 
schemes that enable us to write the above problem in a tractable convex format and finish the section by 
presenting a short procedure that outUnes steps needed to solve the robust WSRM problem. 

A. First Approximation 

In order to arrive at our first approximation of the robust counterpart of the WSRM problem, we will 
handle the uncertainty constraints in (II Ibl i and (II let separately using two different strategies. Let us 
first deal with the constraint in (II Ibl i. Before proceeding we note that for the case of perfect CSI, it is 
possible to ensure the equality constraint on the imaginary part of the desired signal without affecting 
the optimality. However, when the channel is corrupted and the imperfections have to be taken into 
account in the design process as well, the same principle of having the beamformers Wf, ^ orthogonal 
to all channel realizations in the uncertainty set is not a feasible option anymore. Hence, the problem 
is relaxed by dropping this stringent constraint. Nonetheless, it remains to note that owing to the fact 
that Re(c) < |c|, the resultant approximation ensures that if the relaxed problem is solved it also solves 



the original problem. Therefore, by dropping this constraint we obtain a lower bound to the original 
problem. A similar approximation has also been used in the earlier work of lIlTI . Therefore from now on, 
the absolute function in (II Ibl ) is replaced by the real operator that furnishes the real part of the left side 
of the inequality in (II Ibl ). We do not explicitly mention this real operation in the discussion to follow. 

Now to arrive at the tractable representation of the robust counterpart of the inequality constraint in 
(II Ibl ). we need to deal with the following optimization problem 

Pb,k = , min h,,k'^b,k (12) 

where Sb^^k = {^b„k ■ KkK„k^h,k < 1, 9 = 1, • • • , <3} and Z^^^^ ^ P^k'^b^k- the purpose 
of ensuring tractability and equivalence, we rely on a well known result in the duality theory of conic 
optimization. Here it is pertinent to mention that similar approach was also used to obtain a tractable 
formulation of a linear program with polyhedral uncertainty affecting its parameters ||26l . Consider 

mill Re(f"x) : ||Aix||2 < di, Mi (13) 

where f G C", Aj G C"'^",bj G C^^di G M represent data and x G C" is the decision variable. The 
dual of ( [T3] ) can be written as l|23l 

max -A^d : f + ^ Af Ui = 0, ||ui||2 < A^, Vi (14) 

i 

where Aj G M and Uj G C"' for all i are the dual optimization variables. Now if there exists a xq such 
that ||AjXo||2 < di holds (Slater's constraint qualification condition), ([T3] l and (IT4l ) have the same optimal 
values. Using this result, we observe that ([T2] | and the following problem 

max -J]A^_,:w«, = -j;u^X,fc' \Kkh ^ (^5) 

where Z^^ ^, = y'z^^ ^, have the same optimal values when the Slater's condition is valid. Equipped with 
this result, we observe that the equivalent uncertainty immune version of (II Ibl ) can be written as 

Km-Y.Kk>f^''^ii^~^)--<k = -Y.<Ak^ \Kkh<Kk^ vg (16) 

which is clearly a tractable formulation of the inequality constraint in (II Ibb . 

Now that we have dealt with dllbK let us treat the constraint dl lc| ) in the robust counterpart of the 
WSRM problem. For this purpose we note that, after introducing additional variables, the constraint can 



be equivalently written as a set of the following constraints 

B 

j&Ab\k n=l,n^b 



(17) 



\{hb^,k + ^b,,k)'^b,j\'' < I3l^k-,y^b„k ^<Sb„k, \0^b„,k + ^b„,k)'^n,i\^< $bi'^^b„,k £ Sb^^k,n b. (18) 



Clearly, in the above formulation, constraints in (USI l are the troublesome ones. They can be rewritten as 

max \{hb^,k + ^b„k)'Wb,j\'^ < max \(hb^^,k + ^b„,k)^n,lf< l^bt- (19) 

To deal with the left side of inequalities in ([T9t . it is imperative to consider exploring approximations. 
For safe approximations of constraints like ([T9] l. techniques based on Lagrangian relaxations have been 
developed and extensively studied in the optimization literature, see ||T2| and the references therein. 
Borrowing similar ideas, we have outlined a procedure in the Appendix that briefly sketches a proof of 
deriving an approximate LMI representation of the uncertain quadratic constraints under consideration. 

After having outlined the methods needed to arrive at tractable representations of the uncertain con- 
straints given in (lllbl ) and (11 let , we are now in a position to present our first tractable version of the 
robust WSRM problem that can be formulated as 



maximize 



n n " 

\beB keUb 
subject to hbi,^k^b,k 



soc 



bb,k^ \\^b,k\ 



B 



2<\U, yq,beB,keUb 



j&Ab\k n=l,n^bleU„ 



^K:k > 



/5L 



jq b.k 



l^q '^b,k^bb,k 



-^bb,k^b,j 



1 



h'b.k l^q ^b,k 







l^q'^b,k ^b„,k 



E W^bMll < Pb V6 
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where is given in (H, w^.^ G C^,u^,fc G C'^"^ , {tb,k, f^b,kJb,kJli ^1^, ^t^k'^'} G M+ and 

^ G M are the optimization variables. To deal with the constraints in (fT9l ) we have recalled the result 
derived in the Appendix (cf. (02])). Here we note that in the above representation the worst case complexity 
of the above problem will be dominated by LMI constraints [23 |. 

B. Second Approximation 

In this subsection, we attempt to approximate the uncertainty set by an approximately equivalent single 
ellipsoid. Once this is accomplished, a straightforward application of the 5-procedure should reveal robust 
version of the WSRM problem. The problem of approximating complex sets with ellipsoids has been 
studied in various different contexts mainly in the literature pertaining to control theory (see ifTSl for 
references). Ellipsoids are mostly considered for such an approximation because of their simple and 
elegant mathematical properties. It remains to note that these approximations are mainly divided into 
categories of inner and outer ellipsoidal approximations. In both cases, the design problem is to find the 
parametric description of the best ellipsoid that accurately describes the set by either remaining inside 
or outside of it. By the best ellipsoid it could mean to find an ellipsoid of maximum and minimum 
volume in the inner and outer approximations, respectively. Similarly, instead of volume, maximizing the 
minimum axis length can also be used as a criterion for designing the ellipsoid. However, based on the 
structure of the set to be approximated it is not always possible to find both inner and outer ellipsoidal 
approximations. For example, following outer approximation philosophy, a minimal volume ellipsoid can 
be determined to cover a convex hull of ellipsoids but the problem is intractable if inner approximation 
technique is used to determine a maximal volume ellipsoid. For details pertaining to this issue the reader 
is again pointed to the text ifTSl . |[23l and the references therein. 

In order to arrive at the second approximation we first recall a few fundamental results. Let us consider 
the following general representation of an ellipsoid 

£^(E, c) = {x : (x - c)"E(x - c) < 1} (21) 

where E : det(E) > defines the ellipsoid centered at c. The volume of this ellipsoid is proportional to 
(det(E))"^/^. Let us set (E)"^ = E'E'". With u = E'"^(x - c), an equivalent definition of the above 
ellipsoid is 

S'(E', c) = {x = E'u + c : u"u < 1}. (22) 

We note that we can arrive at the same ellipsoidal description as given in (l22l) if we premultiply the matrix 
E'^^ with an unitary matrix U. Based on this non one-to-one behavior the above ellipsoid is referred 



to as 'flat' |[27l . The model in (|2T]) is often more useful when the ellipsoid is not ill-conditioned and 
is nondegenerate ||27| . Now let us deal with a question of prime interest for us. We aim at determining 
the conditions under which the ellipsoidal set of the type defined in (|22)) is contained in the set defined 
in ([2TI ). This is precisely what ensures that £''(A',a) C £ {(DI)^)-^ , d) holds true. It is shown in HH, 
|[23l that the above subset inclusion relation is valid iff there exists A > such that 
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Equipped with the above result let us consider the following set 
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where £i = iS(Ej,Cj) = {x : (x — Cj)^Ej(x — Cj) < 1}. Clearly, Sin represents an intersection of e full 
dimensional ellipsoids of the form given in (1211 . Our problem of interest is to find parametric description 
of a set that would accurately approximate the set <Sj„. To arrive at the optimal parameters of the inner 
approximating ellipsoid, we make the simple argument that if the approximating ellipsoid <S'(E^,Ca) 
is to be a subset of Sin, it implies that <S'(E'„,Ca) C <5(Ej,Cj) for all i. Hence, using volume of the 
approximating ellipsoid as the parameter describing the closeness of the original and the approximating 
set, the approximating ellipsoid can be obtained from the following optimization problem 
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where we have used the fact that for all i, Er^ = DjD^. The above problem is an SDP and it maximizes 
the volume of the inner approximating ellipsoid with E^, as decision variables. The best approximating 
ellipsoid is thus given by {x = E'„*u + Ca* : uu^ < 1}, where E'^j* and c^* are the optimal solutions of 
the above SDpj] Here we recall a result due to Lowner-Fritz John (LFJ) lITSl . ll23l pertaining to extremal 
ellipsoid representations. Once a maximal volume inner approximating ellipsoid of a symmetric set has 
been determined, inflating the ellipsoid by a factor equal to the dimension of the vector space in which it 
is defined {VT in our case), we end up obtaining an outer ellipsoid that contains the original set. The LFJ 



''The case when £i are not inveitible is very similar to the one discussed, and the interested reader is referred to 0231 . 



ellipsoid will be useful in obtaining more conservative approximations of the original robust optimization 
problem. Now armed with the above results, we are ready to describe our second approximation. We will 
devise the new approximation for (II let . The other perturbed constraints in (II Ibb can be handled as in 
the first approximation, albeit with the original uncertainty set replaced with its inner maximal volume 
approximating ellipsoid or the LFJ ellipsoid. Our approximation is build upon the ideas presented above. 
In particular, the approximation is based on the following philosophy 



{5fe„,fc|65„,fcZ^ k^t,k < 1> VQ} « S^^k = {5fe„,fc|66„,fc = uE6„,fc + Cfe„,fe : uu" < 1} (26) 



where ^ = p^^;.?^ and Ej,^^ ^,0;,^ ^ are the parameters of the maximum volume ellipsoid inside 
the intersection of original Q ellipsoids. For a more conservative design, the set ^ may alternatively 
be S^^j^ consisting of a LFJ ellipsoid. The nature of the approximation for Sb^^k is similar. As we detail 
below, the parameters of the approximating ellipsoids can be obtained by solving an SDP similar to the 
one given in (|25] ). We further note that since u = (5;,^ ^ — c^^ fc)(Efe^ fc)^^, we have 

SK,k = {^b„,k ■■ (5fe„,fc - Cb,^,k)'El,kih.,k - Cb^,kf < 1} (27) 

where E^ ^ = (E^ ^Ej,^ fc)~^. Once an approximate description of the uncertainty set like the one 
given in (|27] ) has been established, it is then quite straightforward to obtain tractable version of the 
robust counterpart of the WSRM problem. First we observe that dl lc| ) can be rewritten as 

B 

E M,k+ E ^ ^''',k, \fbeB,keUb (28) 

j&Ub\k n=l,njtbl£U^ 

max |(hfe,_fc + 5fe6,fc)w6j| </3^fc, max \{hb^^k + ^b„,k)'vvn,i\ < $b'k- (^9) 

Here we remark that the above equivalence step is the same as that used in (fT7])-([T8]) and to avoid 
introducing new slack variables we keep the same notation. This should not cause any ambiguity as 
the present approach is independent of the first approximation. However, the uncertainty sets have been 
replaced with the approximate ones that consist of only one best inner approximating ellipsoid. To deal 
with (|29l ). we recall a classical result dubbed as the 5-lemma in the control theory literature ifTSl . 

The simple looking 5-lemma |[T2l . ifTSl can now help us render tractability to the intractable nonconvex 
constraints in (|29] l. Let us without loss of generality focus on the second inequality constraint in (|29] l. 
Clearly, it can be equivalently written as 



where W„ ; = w„ ^wj^^ is the outer product of the vector w„ implying a unit rank constraint on this 
matrix variable. Now a straightforward application of the mentioned above 5-lemma reveals that for 
A6„,fc > the above impUcation is further equivalent to 

-w„,, -w„,,hH, v,^ j-^i. ^..A. y (31, 

Likewise, the first constraint in (|29] l can also be dealt with to arrive at an equivalent tractable formulation 
for the approximate uncertainty set S^j^. Now we are in a position to explicitly state another tractable 
approximation of the robust WSRM problem as 

maximize tf, ^ (32a) 

\beBk(^Ub / SOC 

subject to f dlObl) . (l20cl) for S^^ ^ or S^f,/).^2M. (pOgll, Ab,,fc, A;,,.,^ > : 



eUb\k (32b) 

constraint in (|3T]) \/n G B \ b. (32c) 

Compared to the first approximation, the only new variables introduced are Xb^,k, Ab„,fc G Wj fc, W„ / 
belonging to the cone of positive semidefinite matrices, and we have dropped the rank constraints on 
the matrix variables to ensure tractability of this formulation. We remind the reader that the notation 
(■)soc expresses that the objective can be expressed as a system of SOC constraints. Here we stress an 
important note that dropping the rank constraints, indeed, is a relaxation. However, it has been shown that 
at least in certain scenarios this relaxation is tight ll28l . i.e., once the problem has been solved the rank 
of the matrices obtained is numerically one. Nonetheless, owing to the problem setup in general cellular 
conditions, the relaxation is not tight for the whole range of parameters of interest. Indeed, the recent 
work of |[29l characterizes a similar range. Hence, the beamforming vectors can be extracted from these 
matrices by an appropriate spectral decomposition theorem at least in the case when we have unit rank 
precoding matrices. Now we present the following procedure that can be adopted to yield approximate 
weighted sum rate maximizing beamforming vectors in the presence of channel uncertainties: 



Initialization: set n := and randomly generate 
repeat 

• Solve either the optimization problem in (|20] i for the first approach or (|32] i for 
the second approach (for both cases of simple inner approximation and the LFJ 
ellipsoid based representation). 

• Denote the resulting optimal values of {th^k, fJ-b,k) by (tlk^f^tk)- 

• Set = "P'^^t^ n := n + 1. 
until convergence 



Before concluding this section, we note that it can be mathematically shown that the feasible set of the 
interference terms (II let in the second approximation is a subset of those of the first approximation for 
the case of inner extremal ellipsoid. However, we skip the details due to space limitation and defer the 
comparison discussion to the numerical results section. 

C. Convergence 

We note that the convergence arguments of the two robust approaches and the nonrobust solution are 
very similar. Hence, without loss of generality, we will focus on the first robust approximation presented 
in Sec. IIV-AI We also note that the convergence proof is very similar to the one presented in |[T5l . |l24l . 
Let us define the following set 

CSn = {set of all decision variables in (l20]l|the constraints in (l20bl )-( 20g I are satisfied} (33) 



in the nth iteration of the algorithm that solves (|20] i. Further, let Wn and /{VVn) denote the sequence 
of variables and the objective produced during the nth iteration of the algorithm. In order to conclude that 
f{VVn) < /(^Vn+i), we need to infer some additional intermediate observations. It is clear that Wn 
belongs to both CSn and CSn+i- The inclusion of Wn in CSn is obvious. The inclusion in the feasible 
set of the n + 1st iteration comes from the fact that = ^(t^"'''- 1)/"^ which holds 

because of the conditions given in ^ following the update of variables as mentioned in the algorithms 
for the robust and nonrobust optimization schemes. This in turn amounts to the fact that Wn is contained 
in CSn+i, thus validating our claim. Equipped with this we see that the optimal objective value in the 
n + 1st iteration f{T>Vn+i) is no worse than its value for the variables in the previous iteration i.e., 
/(VVn) < /(^V„+i), hence ensuring monotonicity. Further, owing to the power constraints the cost 
sequence generated in the algorithm is bounded above. Therefore, the proposed iterative procedure is 
guaranteed to converge. The next question of interest is to establish that the point of convergence also 



satisfies tlie Karush-Kuhn-Tucker (KKT) conditions. Since tlie proof strategy for this result is similar to 
the one given in ll24ll . the reader may consult this reference. 



V. Numerical Results 

For illustration purposes, we will focus on a system, sketched in Fig. [T] composed of two cells with 
each cell serving two users. Unless otherwise specified, the number of antennas at each BS is taken as 
T = 4. We will assume that the channel estimates h^^ ^ for all b, k known at the base stations follow the 
CAA(0, 1) distribution. Further, when considering the variation of a quantity of interest with the transmit 
power, we will normalize the transmit power with respect to the noise variance and use the quantity SNR 
instead. Unless otherwise mentioned, we take the weights = 1, for all i. For the parameters considered 
in the simulations, we obtained ranks of the precoding matrices numerically close to one for the second 
approach. This is in agreement with the results presented in |[29l where conditions have been derived 
under which exact rank recovery is possible. Our algorithms usually converged and stabilized within 10 
iterations. 

On the tightness of proposed solutions: It is of interest to compute the gap between the exact results and 
the proposed approximations. Indeed, it appears difficult to characterize this gap analytically. Therefore, 
we resort to numerical tools. To recall, we observe that the exact robust counterpart is given in (ITTi . 
and its approximations are presented in (|20) and (|32] i. To furnish tightest possible approximation to 
([TTI ) numerically, we will need to approximate the exact uncertainty sets with their discrete counterparts. 
Such discrete uncertainty sets will then be used in conjunction with ([TOb to obtain robust solution of the 
WSRM problem based on the worst case strategy. It is easy to see that the sampled versions of the exact 
uncertainty sets would be the subsets of the exact ones, and, hence, the optimal solution would form a 
bound on the exact theoretical solutionis Therefore, we consider a sampled (discrete) counterpart of (ITTI) . 
where the constraints in (11 Ibj i and (11 let are satisfied only for a finite number of samples of the channel 
errors sampled from the uncertainty sets. For mathematical convenience, corresponding to the constraints 
in dllcl ). let Sb^^k be a sampled set of 5f,^ fc, i.e., each element in ^ is also a member of 5;,^ fc- A 
similar notation can be devised for the constraints in (11 lb) . Clearly, the quality of the proposed solution 
would rely on what type of channel errors instances are sampled from the uncertainty sets. Fortunately, 
the tool in ||30l is especially devised for the efficient sampling of the uncertainty sets, and is employed 

'"We remark that on account of finite precision and memory of computers, it is not possible to find tlie exact numerical 
solution as the original uncertainty sets are continuous in nature. 



for this purpose in the paper. To further ensure that the proposed sampling results in tight approximation 
to the exact worst case robust counterpart, we repeat solving the discrete version of (Hi) for a number 
of randomly generated uncertainty sets Sh^,k and Sb^^k- The resulting objective corresponding to the ith 
such run is denoted t*. The constraints in the original problem in (ITTi also satisfy the constraints in its 
sampled version, t* < t* for all i, where t* is the solution of the exact robust WSRM problemllJ Thus, 
the empirical worst case sum rate of the robust counterpart is taken as the minimum of the sequence of 
objectives returned in each iteration. 

In the first numerical experiment, for a given set of estimated channels, the channel uncertainty set for 
each channel from the base station n to user {b, k) is taken as a box of dimension y^, i.e., | [6b„,fc]i| < ^/p 
for i = 1, 2, ... T. To obtain the empirical worst case sum rate, we solve the sampled version of (ITTI ) 
50 times, each with = 10^ independent samples for all n and (6, k). In Fig. |2l we plot the worst 

case sum rate of all the schemes as a function of p. As can be seen, the second proposed design offers 
higher worst case sum rate, compared to the first design, for all values of p in consideration. However, 
once LFJ ellipsoid is incorporated in the second approach, as expected, a decrease in the worst case sum 
rates is observed owing to the conservative nature of the design. Further, we observe that there is indeed 
a gap between the discrete version (labeled as "Empirically, exact RC" in Fig. [2l), and the proposed two 
approaches. The gap remains almost constant for the range of p considered. Minimizing this gap as much 
as possible constitutes a rich area for future research. We also note that we were able to obtain graph till 
p = 0.25 mainly on account of the fact that higher values of this parameter would need unacceptably large 
simulation time to show the variation of the exact worst case robust rate. Nonetheless, such a variation of 
the worst case rates for larger p is shown for the two approximations in the results to follow. For the sake 
of complete comparison we have also illustrated the worst case sum rate achieved by a zero-forcing type 
scheme adopted from ||3], iBl for the scenario under consideration. As expected, owing to the absence 
of its ability to handle channel uncertainties, the rate returned by this scheme remains constant, albeit 
below the nonrobust scheme, for the range of p considered. 

Average worst-case WSR: In the next simulation, we evaluate the average worst case sum rate of the 
proposed robust designs for the given uncertainty sets. We consider the intersection of 3 randomly 
generated ellipsoids for all channels, i.e., Q = 3 in ([4]) to model the uncertainty set. Fig. |3] plots the 
average sum rate of the proposed robust designs versus p for two types of uncertainty sets i.e., (i) 



"indeed, it is easy to infer this from tlie fact tliat tlie sampled version of dllt needs to satisfy fewer constraints than its exact 
continuous robust counterpart. 



the box uncertainty set with deteministic intersecting ellipsoids, and (ii) the above randomly generated 
uncertainty set produced due to the intersection of full dimensional ellipsoids, referred to as complex 
uncertainty set. Furthermore, we ensure that the complex uncertainty set is contained inside the box. 
Thus, as expected, the worst case sum rates of the proposed approaches for the case of box uncertainty 
set are smaller than the complex uncertainty set, which is clearly seen from Fig. |3] It is seen that the 
LFJ ellipsoid based second approximation is the most conservative for the two types of uncertainty sets 
considered above. In fact, the LFJ approximation ceases to remain feasible beyond a certain p. 

To further explore the effect of number of transmit antennas on the performance of the proposed robust 
approximations, we report the above results in a graph for two sets of antennas as shown in Fig. IH We 
observe that larger number of antennas have a more pronounced effect on the average achievable rates of 
the second approach. In fact, the performance of the first approach with T = 8 antennas is very similar to 
that of the second approach with T = 4 BS antennas. This advantage of the second approach disappears 
when we replace the approximate ellipsoids with their more conservative LFJ based counterparts. 

In Fig. [5j we investigate the average worst sum rates as a function of base stations' transmit power. 
In this simulation setup, we consider only set (ii) above, and p is assumed p = 0.1. We observe that, 
compared with the nonrobust solution, the worst case sum rate does not scale significantly with the SNR. 
A possible reason for this behavior could be the fact that to ensure the sum rate problem's constraint 
are met for all channel realizations in the uncertainty set, the base stations have to puU back in terms of 
actually utilizing power for higher spectral efficiency. The conservative nature of the LFJ ellipsoid based 
second approximation is clearly depicted by the bottom most curve of the figure. 

Robustness of the proposed approaches: For the set of nominal channel estimates used in Fig. |2j and 
uniformly randomly distributed errors in a box of dimension p = 0.2, the cumulative distribution function 
(CDF) of the sum rate is illustrated in Fig. |6]at SNR = 10 dB. To gauge the robustness of the proposed 
approaches we calculate the probability that the true worst case sum rate exceeds the objectives (PE) 
obtained by solving (l20l ). (l32l ) and (fTOl ) for the nonrobust design. From Fig. |6] it is evident that this 
exceedance is almost sure for the two robust solutions, while it is a very small number in the nonrobust 
case. Furthermore, although the LFJ approximation is guaranteed to yield PE = 1, it comes at the price 
of lower achievable rate. We further note that the zero-forcing scheme of lH, ||4l, not only yields a 
lower median rate but also totally fails to exceed the PE. In fact, this behavior of the interference nulling 
scheme renders it worse performance than the nonrobust solution of the WSRM problem. In practice, the 
channel errors need not be inside the proposed uncertainty set. Therefore, we introduce a dummy radius 
p' > p and generate Table H The table calculates the PE obtained by solving (|20l ) and (|32l ) for fixed 



p = 0.02, and for uniformly randomly distributed errors in a box of size p' / p at this p = 0.02. It is seen 
that as p' / p increases the second approach shows a more decrease in the value of PE. In comparison, 
when an LFJ ellipsoid is combined with the second approach, owing to its more conservative nature, we 
are gauranteed to meet the PE threshold. 

Finally, we determine a proper value of p in the worst case robust designs, which can guarantee that 
the worst case sum rate is achieved with a given probability. The channel errors ^b^,k are modeled as 
CA/'(0, cjI) for all pairs of k) and n. In Fig. 7(a)[ we calculate such values of p that ensure PE > 







as a function of a for a rectangular uncertainty set in the first robust design approach. Corresponding 
to these p the resulting sum rates and PEs (indicated on top of the markers) for the robust designs are 
provided in Fig. |7(b)| For example, when a = 0.2, we observe that Fig. 7(a) gives p = 0.016. Now 



with this p as design parameter and box uncertainty set. Fig. |7(b)| shows that the worst case sum rate 
(bps/Hz) and the PE are 4.45, 5.8, 3.3 and 0.85, 0.45, 0.98 for the first approach, inner approximation of 
the second approach and the LFJ ellipsoid based second approximation, respectively. It can be inferred 
from this experiment that for a feasible problem, the second approximation combined with LFJ ellipsoid 
provides the minimum rates with almost 100% PE. While the naive inner approximation yielded by the 
second approach results in the highest spectral efficiency (and the lowest PE), the first approximation 
may be considered as a good compromise between the achievable rates and the PE. 

VI. Conclusion 

In this paper, we have studied WSRM problem in a multicellular system by taking into consideration 
the fact that the channels are not perfectly known to the centraUzed base station. Assuming that the 
uncertainties affect the true channels in an affine manner, we design beamformers that maximize system 
wide rates based on the worst case robust optimization strategy. The problem is intrinsically nonconvex, 
NP-hard and intractable. We elevate a novel sum rate maximizing algorithm in perfect CSI to incorporate 
channel imperfections. To do so, we resort to approximating the exact uncertainty set with a tractable 
subset, and thus arrive at two approximations with a varying degree of robustness. The first approximation 
employs a kind of Lagrangian relaxation scheme to arrive at a tractable formulation. The second approx- 
imation relies on modeling the given uncertainty set with extremal ellipsoids. Finally, in the numerical 
experiments we demonstrate the effectiveness of the proposed approaches for different uncertainty regions 
against channel uncertainties. 



Appendix 

In order to arrive at a tractable version of the robust counterpart of the problem in the first approximation 
we, without loss of generality, need to find a tractable representation of the following uncertain quadratic 
form for all b £ B and k i.e., 

^b„,k < ^^b„,k, (34) 

\K„,kWn,lf< ^b^,k, y^b„,k = ^b„,k + ^br„k-\\^b„,k'Pl^^^k\\2^Pbr,,k, g' = l,---,Q (35) 

where we assume that ^ and ^ :^ 0. The condition J2q k ^ ^ implies that the 

uncertainty set defined above is bounded. Recall that for a bounded set C, there exists a number ^ such that 
the distance of all points in C from the origin is bounded above by ^. Indeed, J2q ^6„,fcPfe ^^1^ k — QPbr^,k 
for all 6 € ^ and k G U},. It is clear from (lllcj l that we need to deal with (|35] ) to derive an uncertainty 
immune version of the WSRM problem. (1351) can be equivalently rewritten as 

|(hb„,fc + 5b,^^fc)wn,i|^< ilfe„,fc, y\\h,,,kPll'^k\\l< Pb^,k, q = ^,--,Q (36) 

{^b„,kFl k^f k ^ PKM^ q = l,...,Q^ 
(37) 
56„,fcW„,,6« + 2Re(hfc„,fcW„,;5« J + hfe„,fcW„,/h^^^^ - nb„,k < 0, 

In (|37] ) we have defined W„_i = w„ /w^^. We further note that if Q = 1, we can straightforwardly apply 

5-lemma and arrive at an equivalent tractable representation of the constraint in the form of an LMI. 

Hence, in this case we need to strive for an approximation. We will first make a noteworthy observation. 

It can be seen that if 5;,^ ^ satisfies (|37] ). then so does — 65^ fc. Hence, (l37l ) can be expressed as 

(3o) 

^K,k^n,i^Z,k + 2t6„,fcRe(hb„,feW„,;5H J + hb^,fcW„,ihH - Qb^^k < 

where Z^^_^ ^ PK,k'PL,k- 

Now consider the following relaxation of (|38T l 

Sb„^k^nXr.,k + ^K,kMh^,k^nXr.,k) < {^b„,k - hb„,k^n,l^t,k " Yl ^<l^^Lk + 

q 

T.^'^^b,.,kZl^k€^, (39) 

where Ag > for all q. It is easy to observe that for the conditions stated in (|38T l. the above inequality 
furnishes the implication in (l38T l. Indeed, we see that 

(f)fe„,fc - hb^,kWnMr.,k - Yl \)^K,k + ^1^b„,kZl^^^t,k < ^b^,k - K,k^lVnM^,k- (40) 

9 9 



Hence, it can be concluded that if a tuple (tb„,fc, 6b„,fc) satisfies ( [39l ). it also satisfies (|38] ). Based on the 
above arguments, we can conclude a desirable fact that an optimal solution of the proposed relaxation 
will be a feasible point of the original worst case robust counterpart of the uncertain constraint. Now we 
resort back to our proof and to proceed ahead, note that (|39] is equivalent to 

3A > • I ^^"'^ ~ ^b^,k^nM^,k - T.q\ -hfe„,fcW„,/ \ ^ ^^^^ 

An application of Schur's complement lemma reveals that (1411 can be cast as 



3Aq > 



\-{hb„,k^n,l)^ ^ll 1 



>= 0. (42) 
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Fig. 1. Illustration of a 2-cell system model with 4 users. The dotted-dashed red lines indicate the inter-cell interference, while 
solid black lines show the broadcast part of the signal transmitted by each BS. 



TABLE I 

PE OF DIFFERENT APPROACHES FOR p = 0.02. 
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Fig. 2. Worst case sum rate of different scliemes as a function of p for box uncertainty set. Tfie zero-forcing strategy is adopted 
from ID, El at SNR = 10 dB. 




Fig. 3. Average worst case sum rate of proposed robust designs for different types of uncertainty sets at SNR = 10 dB. 
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Fig. 4. Average worst case weighted sum rate of proposed robust designs for box uncertainty set with T = 4, 8 and SNR = 10 
dB. The weights are taken as (01,1,01,2,02,1,02,2) = (1-14, 1.21,0.52,0.84). 




6. CDF of the sum rate of three approaches when channel errors (for each Hnk) are uniformly distributed in a box of 
0.2 and SNR = 10 dB. The zero-forcing strategy is adopted from ||3], ||4]. 




(a) p that ensures PE is at least 80% versus the variance of channel errors. 
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(b) Resulting sum rates versus p obtained from Fig. |7(a)| The PEs are shown on top of the marker symbols 
used for both approaches. 



Fig. 7. A comparative study of the performance of robust designs with channel errors following Gaussian distribution. 



